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Abstract 

Consider the difference equation = x/ l e a ~ x "-' 1 where a is a positive parameter and d is a 
nonnegative integer. The case d — was introduced by W.E. Ricker in 1954. For the delayed 
version d > 1 of the equation S. Levin and R. May conjectured in 1976 that local stability of the 
nontrivial equilibrium implies its global stability. Based on rigorous, computer aided calculations 
and analytical tools, we prove the conjecture for d = 1 . 
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1 Introduction 

In 1954, Ricker ET1 introduced the difference equation 

x k+l =x k e a - x " (1.1) 

with a positive parameter a to model the population density of a single species with non-overlapping 
generations. The function R\ : K 3 x >->■ xe a ~ x E K is called the 1-dimensional Ricker map. R\ has 
two fixed points: and a. It is not difficult to show that x = a is stable if and only if < a < 2, and, 
for < a < 2, x = a attracts all points from (0, °°) ; or equivalently, the equilibrium x = a of equation 
(jTTTJ is globally stable provided it is locally stable. 

In 1976, Levin and May lfT4l considered the case when there are explicit time lags in the density 
dependent regulatory mechanisms. This leads to the difference-delay equation of order d + 1 : 

x k+1 =x k e a - Xk -«, (1.2) 
* Corresponding author. Email: krisztin@math.u-szeged.hu 
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where d is a positive integer. 
The map 

R d+l : R d+l 3 (x , . . . ,x d ) T ^{x u ... ,x d ,x d e a - x °) T G R d+1 

is called the (d+ 1) -dimensional Ricker map; here T denotes the transpose of a vector. R d+ \ has 2 
fixed points in R d+1 : (0, . . . , 0) r and (a, . . . , a) T . Levin and May O conjectured in 1976 that local 
stability of the fixed point (a, ... , a) T G R d+l implies its global stability in the sense that all points 
from R d + +1 := (0 ,°°) d+l are attracted by (a, . . .,a) T . As far as we know, the conjecture is still open 
for all integers d > 1. The aim of the present paper is to prove the conjecture for d = 1. 

Levin and May's conjecture and many other numerical and analytical studies suggested the folk 
theorem that 'The local stability of the unique positive equilibrium of a single species model implies 
its global stability' This claim was recently disproven by a counterexample of Jimenez Lopez lfTTl[T2l 
on global attractivity for Clark's equation [4] when the delay is at least 3. 

Liz, Tkachenko and Trofimchuk 1 17| proved that if 

0<a< WTv) (13) 

then the fixed point (a,...,a) T G W l+l of Rd+i is globally asymptotically stable, where globally 
means that the region of attraction of (a, . . . , a) T is IR^ +1 . They also suggested that condition ( |1.3| > 
can be replaced by 

0<a< 2^ + 2T^TTF (L4) 

which was proven by Tkachenko and Trofimchuk in [24]. This result is a strong support of the 
conjecture of Levin and May, and it is proven for a class of maps, not only for R d +\- For the 1- 



dimensional Ricker map R\, condition ( |1.4[ ) with d = gives the region < a < 2. For d = 1, i.e., 
for the 2-dimensional Ricker map Rj, condition ( L4) is equivalent to < a < 0.875. See also ifTol 
and lfT5l in the topic. 

Linearising R2 at the fixed point (a,a) T shows that local exponential stability of (a,a) T holds for 
< a < 1, and (a, a) T is unstable for a > 1. As a passes the value 1, a Neimark-S acker bifurcation 
takes place at the nontrivial fixed point. In this paper we show that global asymptotic stability is true 
also for the parameter values a G [0.875, 1]. We emphasise that our result implies global stability at 
the critical parameter value a = 1 as well. 



In case d = 1 the difference equation ( 1.2 1 is equivalent to the 2-dimensional system 



x k+l =y>k, 

(1.5) 

yk+i =yue a Xk , 

and this is also equivalent to the 2-dimensional discrete dynamical system generated by the 2-dimensi- 
onal Ricker map R2. As d = 1 will be fixed in the remaining part of the paper, we shall use the notation 
F a instead of R2 =Rz(ct). F a has two fixed points (0,0) r and (a, a) T . From now on, we shall analyse 
the dynamics generated by F a in the positive quadrant lq_ = (0,°°) x (0,°°). Note that F a (R+) C 
andF a (R 2 + ) cR\. 
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We shall use a combination of analytic and computational arguments. The latter is done using 
interval arithmetics, that is a standard in the area of validated or rigorous numerics. Instead of calcu- 
lating with numbers, we use intervals to control the errors introduced by the computer. After finishing 
a computation, the information we obtain is that the true result is contained in the result interval. We 
shall draw conclusions from that. An interval [a] is represented as a pair of endpoints [a~ , a + ] . Having 
a set S or a number r, we denote their interval enclosures by [S] and [r] , respectively. The reader is 
referred to Moore H3, Alefeld Q, Tucker (2U, EfH, Nedialkov, Jackson and Corliss |20l for further 
details. 

The structure of the paper is as follows. In Section [2] we construct a compact region S, which 
is a closed square around (a,a) T , having the property that F a (S) C S and every trajectory enters it 
eventually. In Section|3]we construct an attracting neighbourhood of the fixed point (a, a) T . We use 
two different approaches. First, the linear approximation of F a is applied at the fixed point. Naturally, 
the size of this neighbourhood tends to as the parameter a tends to the critical value 1 . Then we 
analyse the normal form of F a used in the study of the Neimark-S acker bifurcation, and obtain a 
uniform neighbourhood in the parameter range a 6 [0.999, 1] belonging to the basin of attraction of 
the nontrivial fixed point (a, a) 7 '. In addition to the standard techniques applied in the Neimark- 
Sacker bifurcation, we need explicit estimations on the sizes of the higher order terms in order to 
get a sufficiently large attracting neighbourhood of the fixed point (a,a) T . In Section |4] we give 
an overview on graph representations of discrete dynamical systems and show how it can be used 
to study qualitative properties of dynamical systems. For possible future applications we formulate 
two approaches for general continuous maps in Euclidean spaces. In particular, the correctness of 
an algorithm is verified in order to enclose non-wandering points. In Section [5] we combine the 
computational techniques of Section [4] and rigorously show that every trajectory of F a starting from 

enters the neighbourhood constructed in Section [3] This proves our main result: 

Theorem 1.1. I/O < a < 1, then (a,a) T is locally stable, and F£(x,y) — > (a,a) T as n — >• 00 far all 

{x,y) T G R\. 

For the sake of completeness, here we give a proof for all a € (0, 1]. The result is new only for 
a € [0.875,1]. 

There is an appendix containing some large formulae used in Section [3] The program codes and 
results of our rigorous computer aided computations can be found on link HI. 

Consideration of a bifurcation of a given dynamical system is usually used to show that some 
phenomenon appears in the global dynamics of the system as a parameter passes a critical value. 
The invention in our method is that we use the normal form of a bifurcation in combination with the 
tools of graph representations of dynamical systems and interval arithmetics to prove the absence of 
a phenomenon for certain parameter values near the critical one. As we want to construct explicitly 
given and computationally useful regions, the key technical difficulty is the estimation of the sizes of 
the higher order (error) terms in the normal forms. We hope that our proof shows that these ideas are 
applicable for a wide range of discrete or continuous dynamical systems, as well. 

Running the program of Denes and Makay [8], which is developed to (nonrigorously) find and 



3 



visualise attractors and basins of discrete dynamical systems, suggests that the conjecture of Levin 
and May stands for the 3-dimensional Ricker model, as well. In order to prove the conjecture in this 
case and also for larger values of d, an additional technical difficulty arises. Namely, first a center 
manifold reduction is necessary, and the construction of an attracting neighbourhood should be done 
on the center manifold. Among others, an explicit estimation of the size of the center manifold will 
play a crucial role as well. 

Notations and definitions 

Throughout the paper some further notations and definitions will be used. N, No, M, C stand for 
the set of positive integers, non-negative integers, real numbers, and complex numbers, respectively. 
The open ball in the Euclidean-norm ||.| and in the maximum norm with radius 8 > around q G 
R" are denoted by B(q;8) and K(q;8), respectively. In Section [3] we shall often use the notation 
Bg = {z G C : \z\< 8} for 8 > 0, where \z\ denotes the absolute value of z G C. For a vector x G C", 
x T denotes the transpose of x. For § = {^i^i) 7 G C 2 and £ = (£i,£ 2 ) r G C 2 let denote the 

scalar product of them defined by (E, , Q = E,\ £i + ^2 £2- Let also a = (a, a) T . 

Let / : @f C W 1 — > W be a continuous map. For k G No, / denotes the &-fold composition of /, 
i.e., f k+l (x) = f(f k (x)), and f°(x) = x. A fixed point p G Stf of / is called locally stable if for every 
£ > there exists 8 > such that \\x — p\\< 8 implies \\f(x) — p\\< E for all n G N. We say that 
the fixed point p attracts the region U C £Fy if for all points u G U, \\f n {u) — p\\— > as n — > 00. The 
fixed point p is globally attracting if it attracts all of Qif and is globally stable if it is locally stable and 
globally attracting. 

2 The dynamics in the first quadrant 

In this section we construct compact squares sj ^ C M 2 ^, i G No around a so that F a (Sj- ) C s\ and 
5-°^ attracts all points of IR 2 for all / G No and a G (0, 1]. Hence an elementary proof of Theorem |l.l| 
is obtained for < a < 0.5. Recall F a (M 2 ) C IR 2 ^. We can illustrate the image (xk+ 1,^+1) of {xk,yk) 
as first going horizontally from (xk,yk) to the diagonal, proceeding upwards if < xt < a, otherwise 
downwards vertically until we reach the value yt+i = yke a ~ Xk . This is shown on Figure[l] 
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Figure 1: Dynamics for x > 0,y > 

Let (xo,yo) £ and < a < 1 be given. Define the sequence (xk,yk)k=o e ^+ by 

(x k+ i,y k+ i) =Fa(x k ,y k ), k £ {0,1,...}. 

Consider the following cases: 

I. < xq < a < y 

Clearly we have a < x\ < y\ and max{xo,yo> ;c i>.yi} < y\. 

II. a < xo < yo 

In this case a <x\ andyi <x\ wifhmax{xo,yoi-Xi>yi} <yo- We distinguish two cases depending 
on yi < a or not. 

III. a < y < xq 

We obtain a < yo = x\ and y\ < yo = x\. During the consequent iterations y,-_|_i < y,- = x i+ i is 
satisfied as long as a < x ; - stays true. If a < x,- for all i, then y; > for all i, and both (yi)°° = o and 
(x,)~ =0 are monotonically decreasing, and converge. The only possibility is (xk,yk) — > (cc,a), 
since the other fixed point is at (0,0). Otherwise there is a minimal N > such that < yw < 
xn < a is satisfied. We note that < yx-i < a < x^_i is true. We have max ie { r „m{xi,yi} < xo. 

IV. < y < a < x Q 

Obviously < yi < x\ < a, and max{xo,yo,xi,yi} < x o- 

V. < yo < xo < a 

Here we have < x\ < a, < x\ < y\ and max{xo,yo,xi,yi} < yi- We distinguish two cases 
depending on a < yi or not. 
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VI. < x < y < a 

Then x\ < a and x\ <y\. Now xt+i = yt < y ; -+i is satisfied as long as Xi < a stays true. If 
Xi < a for all i, then y; > for all i, and both {yi)°°=o and (jc,-)^L are monotonically increasing, 
and converge. The only possibility is (xk,yk) — > (a, a), since the other fixed point is at (0,0). 
Otherwise there is a minimal N >0 such that a < x^ < yN is satisfied. We note that < x^-\ < 
a < y N -i. We have max, G { () .... !W }{^,y,} < y N . 

We conclude that for any (xo,yo) G the sequence ((jCjt,)'jt))r=o converges to (a, a) or enters the 
triangle H{h,Q ) = {(x,y) : < y < x < a}, with = 0, that is the region denoted by V in 
Figure [2] We mark the possible transitions with arrows, if it is solid, it refers to transition in one step, 
otherwise it is possible to have multiple iterations before entering the next region. 




Figure 2: Dynamics for x > 0,y > 

Assume now that (xo,yo) € H(h ). Either the sequence converges to (a, a) or there exists an 
N > I such that (jtjv-i,yjv-i) £ VUVI, (xjv,yjv) £ I an d (xN+i,y>N+i) £ //. So the following stands: 

y N+l =y N e a - XN =y N ^e 2a - x »- XN - i < ae 2 ^ 21 ^ . 

This implies that (x N+l ,y N+i ) £ G(g^), where G(g^) = {(x,y) : a < x < y < g^}, with g^ = 
ae Za ~ 2h o . Now there exists M > N + 2 such that (x M -i,yM-i) 6 II U III, {xm^m) G ^V and 
(im+i ,yM+i) G V. We get the following inequality: 

y M+l =y M _ ie 2a - x M-^-i > ae 2a-2gtf\ 

Therefore (xM+i,yM+i) £ H(h[ a ^), with /i^ = ae 2oi_2s » . Similarly, the sequence will enter the 
set G(g[ a ^), with = ae 2a ~ 2h ^ . Repeating this argument we get two sequences (/z,- a ')°l and 
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{gi a) )?=o defined recursively by h { Q a) = 0, h\ a) = ae 2a " 2 «w for i > 1, and gi = ae 2 "-^'" 1 for z > 0. 
It is easy to see that (/i|° < ' > ) increases and (g\ ) decreases. We have the limits /iL°^ < a and 

> a. Define 

s i ={{x,y)- h i <x<gi A- <y<8\ }• (2.i) 

We sum our observations in the following Lemma: 
Lemma 2.1. For every a G (0, 1] and for every i G No 

2. s[ a ^ attracts all points o/M+. 
(cc) 

If Soo = {(a, a)}, then we have already established that the fixed point is globally attracting. 

(oc) (a) 

This is not the case however for a > 0.5. For a given a, we can view the sequences (h\ ) and {g\ ) 
as even and odd iterates of the function 

T«(f) = ae 2a - 2t , (2.2) 

starting from t = 0. That is, h [ Q a) = T°(0), g ( a) = T a (0), h[ a] = t£(0), g[ a) = ^(0), ... . The non- 
rigorous bifurcation diagram on Figure [3] shows that the unique fixed point a of T a is stable for < 

(oc) 

a < 0.5, it is unstable for a > 0.5, and there is an attracting 2-cycle for a > 0.5. Thus SL = {(a, a)} 

(cc) (cc) 

can be expected only for a G (0,0.5). For a > 0.5 the two points of the 2-cycle give hi, ' and gi, '. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure 3: Bifurcation diagram of z a 



Proposition 2.2. I/O < a < 0.5, then for every (x,y) G we have F£(x,y) — > (a, a) as n — > oo. 

Proof Observe that z' a {t) = -2z a (t), z a (t) > 0, and if t > a and a G (0,0.5], then 2z a (t) = 2a<\ 
is satisfied. Since 



dt 1 



(T«(T a (0) -0 = 8T«(T a (0)T«(0 (2t„(0 - 1) , 



(2.3) 
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therefore z a (z a (t)) — t is a concave function for t > a. The first derivative at a € (0,0.5] is 



dt 



(To(Ta(O)-O 



= 4T (T«(a))T„(a)-l =4a 2 -l < 0. 



(2.4) 



These imply that 4 (T«(T a (f)) — f) < for all ? > a. Thus the only zero of T a (%(?)) -f in f G [«,<») 
is a. This gives us that g» = a which implies that /zL"' = a and consequently sL°^ = {(a, a)} in 
this parameter region. □ 



We assume a G [0.5, 1] in the remaining part of the paper. 



3 Attracting neighbourhood 

Let us consider map F a . In this section we are going to give e(a) > 0, such that 
inf ae [ 5 !j e(a) > and K(a;e(a)) belongs to the basin of attraction of a for a G [0.5,1], that 
is, Fa(xo,yo) — ^ a as « — > oo for all (jco,;yo) r G K(a;e(a)). 

Introducing the new variables u=x — a, v = y — a, map F a can be written in the form 

( " ) h^A(a) ( " ) +/ a («,v), (3.1) 



v / \ v 



where the linear part is 



and the remainder is 



A(o) 



1 

-a 1 



fa(u,v) 





v(e- M -l) + a(e-"-l + M) 



The eigenvalues of A(a) are /112(a) = 1±l ^^ a 1 g C, and the corresponding complex eigenvectors 
are q L2 (a) = (j^^^ = ,1^ G C 2 , respectively for a > i. Let 9 (a) = ?i(a) 

and /i (a) = /^i (a). Let jc(a) G C 2 denote the eigenvector of A(a) T corresponding to ju(a) such that 
(p(a),q(a)} = 1,. This leads to 

/ in a /zL/t — 1 -I- j\ 

p(o) = 



/a \/4a - 1 + / 
V4a-l' 2V4a-l 



We introduce a complex variable 



z = z(u,v,a) = ( P (a),(u,v) T ) = \ (v- ^^^i ) • (3 " 2) 



We also have an explicit formula for (u,v) T in terms of z, which reads as 

r 



, ,t , , — 7 -^- l Rez + V4a - llmz \ ,„ „ 

(w,v) r = 24(a) + 29(a) = ( ^ ,2Re z j . (3.3) 



Our original system ( |3- 1[ > is now transformed into the complex system 

z>-> G(z,z,a) = (p(a),A(a)(zq{a)+zq{a)) + f a (zq(a) +zq(a)) 
= H(a)z+g(z,z,a), 



(3.4) 
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where g is a complex valued smooth function of z, z and a, defined by ( |A.1| ). It can also be seen that 
for fixed a, g is an analytic function of z and z and the Taylor expansion of g with respect to z and z 
has only quadratic and higher order terms. That is, 

1 



> z >°0 = L kw 8kl ^ zkzl ' ™& k > l = > l ,-—> 



2<k+I 



where g u (a) 



d k+l 

dz k dz' 



{z,z,cc J 



z=0 



for& + />2, k,l = {0,\,...}. 



Proposition 3.1. Let a E [0.5, 1) be fixed and 
e(a) = min 



1 /4a -1 9(4a-l)(l-Va) 



20 V 2 + a '20(l + 2Va)\/2 + a J ' 



77ze?i j(.x,y) r E M 2 : |jc — a|< e(cc),\y — Ct\< e(oc)} belongs to the basin of attraction of the fixed 
point a ofF a . 



Proof. Let us study the map in the form ( |3.4| >. First note that ( |3.3[ ) easily implies inequalities 

2 

|w|< — and |v|< 2\z\. 
Ja 



(3.5) 



Assuming |w|< 1/10 and |v|< 1/10 and applying the inequalities \e 11 — 1| < e 1 / 10 !^! < ^\u\ and 
\e~ u — 1 + u\ < e l / l °Y < |m 2 we obtain the following estimations: 



\g(z,z,a) 



p(a),f a [zq{a)+zq(a 
sJAa - 1 + i 



v{e~ u -l) + a(e~ u -l + u) 



2V4a-l 

s^(Mk--i|-H»k--i+^)s^ 

5 



|«vk 1/10 + a- 



,1/10 



< 



9 V 4a - 1 v 
Now, since |/i(a)|= A/a, hence 



*+2H))<^- 4a+V5) tf. 



9 A /a(4a-l) 



lrt _ ,7 r 5 4(l+2v/a . A . . . . 
G(z,z,a) < \/a+Q-^==N z < N 
V 9 V«4a-1 y 



provided that IzM is so small that \u\< 4* and |v|< jk and + § • 4 ( 1+2 v / «) : | z | < 1 Inequalities 

' 1 y A/a(4a— 1) 



3.5| ) imply that |z|< ^ guarantees |«|< ^ and |v|< jq. Therefore 

„ , , , a . fv/a 9(l-v/a)\/a(a-l) 



(3.6) 



implies |G(z,z, a)|< |z|, which means that \G n (zo,zo, oc) | — >■ as « — >• oo if |z |< £g(&) is satisfied. We 
show this by way of contradiction. Assume that |zo|< £g(o0> in = G n (zo,z~o,Cc) and \zo\> \zi\> ■ ■ ■ > 
\z n \> • • • > with \z n \— > c > as n — > oo. Since G is continuous we have that max| z | =c |G(z,z, a)|< c 
and consequently |zjt|< c also holds if k is large enough which is a contradiction. 
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From equation (|3.2|) one obtains that if |m| < 5, |v|< 8, then 



From p.7| > we infer that if |h|< e(a) and |v|< e(a) then |z(w,v,a)|< £g(o5) which completes our 
proof. □ 

Note that e(a) goes to as a goes to 1. This means that when a is close to 1, then the constructed 
region K(a;e(a)) becomes very small. Thus it is impossible to show by interval arithmetic tools that 
every trajectory enters into it eventually. Nevertheless s(a) might be used in the case a G [0.5,0.999]. 
However, our following method is not only capable to give an attracting neighbourhood for all a G 
[0.999, 1], whose size is independent of the concrete value of the parameter, but it also serves a better 
estimation of the attracting region even if we assume only a G [0.875, 1]. The main results of the 



section are the following two propositions. Here, we only present the proof of Proposition 3.3 The 
whole argument can be repeated to get a universal attracting neighbourhood when only [0.875, 1] is 
assumed. The differences only appear in concrete values in the given estimations. 

Proposition 3.2. Forallfixeda G [0.875,1], {(x,y) T el 2 : |x-a|< 1/37, |y— Cf|< 1/37} belongs 
to the basin of attraction of the fixed point a ofF a . 

Proposition 3.3. Forallfixeda G [0.999,1], {(x,y) T el 2 : |x-a|< l/22,\y-a\< 1/22} belongs 
to the basin of attraction of the fixed point a ofF a . 

Proof. We follow the steps of finding the normal form of the Neimark-Sacker bifurcation, according 
to Kuznetsov lfl3l . In our calculations and estimations we use symbolic calculations and built in 
symbolic interval arithmetic tools of Wolfram Mathematica v. 7 or 8. 

According to Kuznetsov lfl3l . we are aiming to transform system \2>A\ to a system which takes 
the following form. 

w i y iu.(a)w + c l (a)w 2 w + R 2 (w,w,a), (3.8) 

where c\ and R 2 are smooth, real functions such that for fixed a, R 2 (w, w, a) = 0(\w\ ). We are going 
to show that there exists So > such that for all | w| < £o and a G [0.999, 1], 

|ju(a)w + ci(a)w 2 w + R 2 (w,w,a)\ < \w\ 

holds, which implies that B £o belongs to the basin of attraction of the fixed point of the discrete 
dynamical system generated by ( |3.8| ). From this, we shall be able to show that the fixed point a of F a 
attracts all points of K(a; ^). 

To be more precise, for a fixed a, we are looking for a function h a : C — > C, which is invertible 
in a neighbourhood of G C and which is such that in the new coordinate w = h^, l (z), our map ( |3.4| ) 
takes the form 



w^-h a l (G(h a (w),h a (w),a)) = ii(a)w + ci(a)w 2 w+R 2 (w,w,a), (3.9) 
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where R 2 {w, w, a) = 0(|w| 4 ) for fixed a. According to |[T3l , such a function can be found in the form 

/120(a) 2 _ /102(a) _ 2 , %)(a) 3 , /212(a) _ 2 , /io3(a)_ 3 ^ im 

h a (w) = wH ^—^-w +hn(oc)ww-\ ^ — -w H -f-^-w -| ww H -^w J . (3.10) 

2 2 6 2 6 

Clearly, /i« has an inverse in a small neighbourhood of G C. A formula for /i^ 1 can be given in 
the form 

fcaU-z) =h im , ;a (z)+R3(z,a), 

where 

2<k+l<4 

and 7?3(z,a) = 0(|z| 5 ). The coefficients can be obtained by substituting w = h~ l {z) into z = h a (w) 
and equating the coefficients of the same type of terms up to the fourth order. The result for hi nv ,a 



in terms of /220(a), . . . , /203(a) is given in ( |A.14| ). The coefficients /120(a), . . . , /103(a) are determined 
such that 

h a l (G(h a (w),h a (w),a)) 

has the form ju(cu) + c\ (a)w 2 w plus at least fourth order terms in w, that is, the transformation kills 
all second and third order terms with one exception. This requires the condition 

ft \ \ k 

) \ ^ t forallAcG {1,2,3,4}. 



As /J. (1 ) = i + i^, this is clearly satisfied. Formulae ( |A.15| )-( [A.20| ) contain the obtained results. 

In which region is the transformation valid? We are going to show that h a is injective on B1/3 C C 
and that /i^ 1 is defined on B^. Let us suppose that z G C is fixed and /120(a), . . . , /103(a) are given 
for a fixed a G [0.999, 1]. Let 

H az :C3tt'4w + z- h a {w) G C. 

By this notation, H az {w) = w holds if and only if h a (w) = z. Now, we have the following 

\H a , z {wi) -Ha tZ (w 2 )\ =\wi -h a (w\) -w 2 + h a (w 2 )\ 
<\w\ — w 2l , . 



■((^ + IM«)I+^)(M+M) + 

(!^ + ^M + ^l) (Kf+MM+M 2 )) 



If |/i2o(a)|/2 + |/iii(a)|+|/io 2 (a)|/2< 5i, |A 30 (a)|/6+ \h l2 (a)\/2 + |/i 03 (a)|/6 < 8 2 , \w\< 83 and 
\z\ < 84 hold, then we have 

\H a ,z(wi) -H a , z (w 2 )\ < \w\ — W2|(25i53 + 3<5i<5 3 2 ) 

and 

\H a , z (w)\ <8 4 + 8 l 8^ + 8 2 8i. 

By interval arithmetics we obtain that the first two inequalities are fulfilled if 81 = 0.76 and 8 2 = 0.52. 
Now, if we choose 83 = ^ and £4 = | we obtain that H az : B1/3 — > B1/3 is a contraction. Hence 
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for all fixed z G B1/5 there exists exactly one w = w(z) G Bj/ 3 such that H a ,(w{z)) = w(z), that is 
h a {w{z)) = z. This means that h a l can be defined on Bi/ 5 . 

The obtained estimations on h a are going to be used in the sequel. These were 

|/j 2 o(a)l/2 + IMa)l+IM«)l/2<0.76, 

|Ma)l/6 + IM«)l/2 + IMa)l/6<0.52 

and in particular 

|w|-0.76|w| 2 -0.52|w| 3 < \h a (w)\< |w|+0.76|w| 2 +0.52|w| 3 (3.11) 
for all a G [0.999,1]. 

Let H = {(a,w) elxC: a G [0.999, 1], w / and \w\< 1 /20}. In the sequel, we shall 
always assume that (a,w) G H. From this assumption and inequality ( |3.11| > we readily get that 
|w|< 1.05|/j a (w)| from which we get in particular that \z\= \h a (w)\< 1/19. 

Our goal now is to give £0 G (0, 1/20], independent of a such that for every a G [0.999, 1], if 
< \w\< £0, then 



h a l (g (h a (w),h a (w),a 



< \w\ holds which guarantees that B eo belongs to the 
basin of attraction of the fixed point of the discrete dynamical system generated by ( |3.9| >. For, we 
turn our attention to the estimation of function R2 in (|3.9|). 



First, we go back to ( |3.4[ ). Let us consider 

" k\l\ 



k+l=23 



The explicit formulae for g%Q{oc), ■ ■ ■ ,g(&{oc) can be found in equations ( |A.2| )-( [AT8] ). By interval 
arithmetics, one may obtain that for all a G [0.999, 1], 

^ + kn(«)l + l^ = i^^^<1.01, (3.12, 
2 2 A /a(-l+4a) 



|g30(«)| |g2l(«)| |g!2(«)| jgtoCgOJ 



(6 + a) ; /2 + (a-2)a 
9a(4a-l) + \/ a 2 (4a-l) 



< 1.09. (3.13) 



We also have that flj (z,z, a) = £^+^4 8j j£ L z k z l +R\ (z,z, a), where /?i (z,z, a) = 0(|z| 5 ) for fixed 



a. For explicit formulae of the fourth order coefficients see equations ( |A.9[ )-( |AT3j ) in the appendix. 

(-u) k \ 



It is clear from equations (|3.1[), (|3.3|> and (|3.4[) that 



/?i(z,z,a) 



V4a - 1 - / 
2V4a - 1 



f (-«*)* v 



/t=4 



/t=5 



where w and v are defined by equation p.3| >. Using < \z\< 1/19 and ( |3.5| ) we have that |m|< 
1/8, |v|< 1/8 and obtain 



\Ri(z,z,ot) \ < 



a 



4a -1 



e i/s 
IT 



,1/8 



< 



a 8 



r 



5! 
4 \ 



\u\ < 



a 8 / . . 16 32 



4a- 17 



24a 2 a 3 /2i20' 



64 



4a - 1 7 V57a 2 + 285a 2 J ' Z ' 665 V (4a - l)a 3 ' 
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Now for all (a,w) G H with z = h a (w) we get that 



\Ri(z,z,a)\< £ 



k+l=A 



gki(u) k-l 
z z 



kill 



+ \Ri(z,z,a) 



6-3a + v / a(12 + a)+4V3 + a 2 
12 v / a 3 (-l+4a) 



+ 7?i(z,z,a) 



< 



4758 - 1995a + 665 ^a(12 + a) + 2660\/3 + a I 



7980 v / a 3 (-l+4a) 
By interval arithmetics we obtain for all a G [0.999, 1] that 

\Rt(z,z,a)\ <0.76|z| 4 . 

We also need a similar estimation on \h~ l (z) \ . Let us recall that z = h a (w) and 

w = h a \z) =h im ,, a (z)+R3(z,a). 



(3.14) 



As hi mM {z) is a polynomial of z and z of degree four (see formulae ( |A.14| )-( lA.20[ )), we denote the 
coefficient corresponding to z k z l by h k j w (a). Calculating these coefficients and using interval arith- 
metics, we obtain that for all a G [0.999, 1] , 



£ |/4l,| < 0-76, £ \hg,\ < 1.06 and £ \h^\ < 1.39. 

fc+/=2 *;+/=3 /t+/=4 



(3.15) 



See formulae ( |A.21[ ) - ( |A.27| ). Let us recall that for all (a,w) G H, \w\< l.05\h a (w)\ holds. Now, for 
the fifth and higher order terms in /?3, first we give an estimation of type \R$ (h a (w) , cc)\< K^\w\ 4 and 
then we get that \R 3 (z,a)\< K 3 l.05 4 \z\ 4 , with z = h a {w). 

From the definition of h im , a , it follows that R-s(h a (w) , Oc) = w — hi m ^ a (h a (w)) is a polynomial of 
w and w and it has only fifth and higher order terms. Let r^ 1 (a) denote the coefficient of R 3 (h a (w) , a) 
corresponding to w k w l . We use a bit rougher estimation for \Y,5<k+i r t^ 1 i 01 -)^^^- Namely, first we 
give the estimations 



|/j 20 (a)|< 1.01; |An(a)|< 0.001; |/j 02 (a)|< 0.51 
|/j 30 (a)|<0.89; \h n {a)\< 0.45; |/j 03 (a)|< 0.89 



(3.16) 



for all a G [0.999, 1]. Now, in Ri,(h a (w),a) we replace w and w by |w|, h nm (a) by the estimates 
given in ( |3.16| ) (for 2 < n + m < 3, (n,m) ^ (2, 1)), and then we turn every — sign into + to get areal 
polynomial 7?3(|w|), with nonnegative coefficients r\ (independent of a) corresponding to |w|*\ If we 
use < | w\ < 1/20, then we get that 

k-4 



£ \ri'(a)\\w\ 

5<k+l<l2 



k+l < £ %\M k < £ r\\ 



w 



5<k<\2 



5<k<l2 



1 

20 



< 1.02|w| 4 . 



This implies that 



\R 3 (z,a)\< 1.05 4 -1.02|z| 4 < 1.24|z| 4 
for all (cc,w) G H, where z = h a {w). It is now clear that 

h inv {z) < |z|+0.76|z| 2 +1.06|z| 3 +2.63|z| 4 

holds for all (cc,w) G H with z = h a (w). 



(3.17) 



(3.18) 



(3.19) 
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Now, we are able to give an estimation on R 2 in ( |3.9| >. First, according to our previous estimations, 
let us define the following real polynomials. 

h max (s)=s + 0J6s 2 + 0.52s 3 , 

G max {s) =s+1.0l/ + 1.09? + 0.76/, 
h^(s) = s + 0.76/ + 1 .06/ + 2.63/ 

and 

48 

V „ f „\ untax _ /~imax _ unax t „\ 
Q[ s ) = L<lk\s) = h im oG oh (s). 

k=\ 

By now, it is obvious that for (a,w) G H, \R 2 (w,w,a)\ < EibU9ifc(l w l)l w l 4 (5g)* 4 holds, which leads 
to \R2(w,w, cc)\ <23.9|w| 4 . For our purposes this approach is too rough - the obtained neighbourhood 
would be too small (?» K(a; 1 /80)) and we could not prove that every trajectory enters it eventually. 
Hence we have to be as sharp as we can in our estimations to obtain as large neighbourhood as 
possible. So, instead of only estimating these functions separately, let us consider the composite 
function h~ l o G a o h a , where G a (z) denotes G(z,z, cc). Now, we are only interested in the fourth and 
higher order terms. Since h a is a known function and we also know functions G a and h~ l up to fourth 
order terms, hence we are able to compute the fourth order coefficients of h^ 1 ° G a o h a , denoted by 

Jcl 

'2 



r^{a), where k+l = 4. By interval arithmetics we show that 



£ |/f(o)|< 1.02. (3.20) 



k+l=4 



See equations ( |A.28| > - ( |A.32[ ) for the formulae. Using this we infer that 

48 / 1 \ k-4 



48 / j \ K-4 

R 2 (w,w,a)\< 1.02|w| 4 +£^|w| 4 ( — <4.6|w| 4 . (3.21) 

k=5 V 20 / 



for all (a,w) G H. 



Now, we turn our attention to c\ (a) in (3.9 1. The formula for ci(cu) can be found in (A.33 1. 



According to [ 13 ] and using inequality (|3.21|) we get the following 



\n(cc)w + c\(a)w 2 w + R2(w,w,a)\ < \w\\n(a) +c\(a)\w\ 2 \ + \R 2 (w,w,a) 
= |w|||/i(a)|+d(a)|w| 2 | + \R 2 (w,w,a)\ < \w\ U/a + d(a)\w\ 2 \ +4.6|w| 4 , 

for all (a,w) G H, where d(a) = ^j^ci(a). Let 

21 / rz, 1 



(3.22) 



R/[{w,a) = |v« + d(a)\w\ - (\fa + a(a)\w 

where a(a) denotes the real part of d(a). 

In the following we are going to prove that |/?4(w, a)| < 0. 1 |w| 3 holds for all (ot,w) G H. First of 
all, the formula for function a is the following 

, , 4 + a(-10 + a + a 2 ) 

a(a) = — ^ (3.23) 

v ' 4a 3 /2 -i + «4 + a 
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It can be readily shown that 



1 < a(a) < 



1 



(3.24) 



holds for all a G [0.999, 1]. Using the definition of d and a, the estimation above and the assumption 
(a,w) G H we get the following. 



a+d(a)\w\ | - (y/a+a(a)\w\ )| = 
(k/(a)| 2 -H«)) 2 )|w| 



a + 2 v / aa(a)|w| 2 +|d(a)| 2 |>v| 4 - (y / a + a(a)|w| 



y / a + 2 v / a<3(a)|w| 2 +|t/(a)| 2 |w| 4 + *fa + a(a)\w\ 2 
(|J(a)| 2 -(«(a)) 2 )|w| 4 



< 



^400a|w| 2 +2Vaa(a)|w| 2 + 20Va|w|+a(a)|w| V400a - 2 + 20a - 1 



< 



(|^(«)| 2 -(a(«)) 2 ) 3 ^ (|J(«)| 2 -(a(«)) 2 ) 3 
= • w — • \w\ - 



{\d{a)\ 2 -{a{a)f 



w 



'a 4 + 3a 3 - 12a 2 + 20a -4) 



19^0 + 18a 



37 a 



16 • 37a 6 (4a - l)(a 2 +4a - l) 2 

From this last formula it can be proven that |/?4(w, a)|< 0.1 |w| . From inequalities ( |3.22| ), ( |3.24[ ) and 
from the above estimate we obtain that for all (a,w) G H the inequalities 



\jj.(a)w + ci(a)w 2 w+R2(w,w,a)\ < \w\ (y/a + a(a)\w\ 2 +R4(w,a)) +4.6|w|' 



1 



< V«k|--|w| +4.7|w| < \w\ 1 - -\w\ (1 -4-4.7|w|) < \w 



(3.25) 



hold provided that |w|< ^ = min { ^ , } . This means that for all a G [0.999,1] the set Bj/20 
belongs to the basin of attraction of the discrete dynamical system generated by ( |3.9| ). 

Now, we only have to show that for any (x,y) T G K(a; 1/22), after our transformations, \w\< 
£ o = 25 h^s. First, we need e G such that for \z\< £g, | w| = \h~ l (z)\< ^ holds for all a G [0.999, 1]. 
From ( |3.11| ) we obtain that Eq = ^ is an appropriate choice. Now, if |«|< e, |v|< e, then from ( 3.7 1 
we get that 



\z(u,v,a)\< e 



Thus for 



£ < 



a(2 + a) 
4a -1 



4a- 1 1 



mm 



ae[o.999,i] V a(2 + a) 21 



(3.26) 



we get that if \u\< e and |v|< e then |z(w,v,e)|< ^j. It is easily shown that e = X= fulfils inequality 



p.26[ ) which proves our proposition. □ 

Proposition 3.4. Fixed point a of map F a is locally asymptotically stable if and only if a G (0, 1]. 

Proof. By linearisation one readily gets that a is locally asymptotically stable if a G (0,1), and 
unstable if a G (1,°°)- Transforming map F a to the form ( |3.9| > and using inequality p.25[ ) yields that 
I is a locally asymptotically stable fixed point of F\ . □ 
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4 Graph representations 



Covers and graph representations 

Different directed graphs can be associated with a given map. The graphs reflect the behaviour of 
the map up to a given resolution. The vertices of these graphs are sets and the edges correspond to 
transitions between them. We can derive properties of our dynamical system through the study of the 
graphs. These techniques appeared in many articles, in both rigorous and non-rigorous computations 
for maps by Hohmann, Dellnitz, Junge, Rumpf J6), 0, Galias [9], Luzzatto and Pilarczyk [QjQ, and 
computations for the time evolution of a continuous system with a given timestep by Wilczak ETTl . 
We introduce the general setting and two applications in particular. One to enclose the non-wandering 
points and the other one to estimate the basin of attraction. 

Both methods (Algorithms [T] and [2]) combined with local estimations of the type of Section [3] at 
the critical points, can be theoretically applied to prove different dynamical properties. On the one 
hand these algorithms are included for possible future references, on the other hand certain elements 
of these algorithms proved to be useful for the map F a in Section [5] In particular, the correctness of 
Algorithm [T] is crucial in Section [5J 

Definition 4.1. .5" is called a cover of S> C W 1 if it is a collection of subsets of W such that U se ys D 
We denote their union U se ys by \y \ in the following. We define the diameter or outer resolution 
of the cover by 

= diam(^) = supdiam(s). 

with 

diam(^) = sup||x — y\\. 

x,y£s 

A cover S? 2 is said to be finer than the cover 5f\ if 

Oi G &i ) (3{s 2 j,i G J?} C y 2 ) such that \J s 2 j = s v 

We denote this relation by S? 2 - The inner resolution of a cover 5? is the following: 

M-{Sf) = sup{r >0:Vx£@,3s£ r y : B(x;r) C s}. 

A cover 5? is essential if 5? \ s is not a cover anymore for any s G .5f . The cover 8? is called a 
partition if it consists of closed sets such that \ 8?\= Q) and ~ip\,p 2 G 8? : p\ Pip 2 C bd(pi) Ubd(p 2 ), 
where bd(p) is the boundary of the set p. Consequently, for any partition & the inner resolution 
<%T{&)S& zero. 

In the following we will always work with essential and finite partitions, as a consequence the 
supremum in the definition of the diameter ^ + of the partition becomes a maximum. 

Definition 4.2. A directed graph <S = §f , <?) is a pair of sets representing the vertices "V and the 
edges S, that is: $ C y x "V , and (u,v) G <§ means that <£ has a directed edge going from u to v. We 
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say that vi — > V2 — > ■ ■ ■ — > Vk is a directed path if (v,-, v,-+i) G # for all / = 1, ... ,A: — 1. If v# = vi, then 
it is a directed cycle. 

A directed graph is called strongly connected if for any u,v£f,v/u there is a directed path 
from u to v and from v to w as well. The strongly connected components (SCC) of a directed graph 
are its maximal strongly connected subgraphs. It is easy to see that u and v are in the same SCC 
if and only if there is a directed cycle going through both u and v. Every directed graph <£ , can be 
decomposed into the union of strongly connected components and directed paths between them. If we 
contract each SCC to a new vertex, we obtain a directed acyclic graph, that is called the condensation 
of#. 

We say that the directed paths p\,p2 are from the same family of directed paths, if they visit the 
same vertices in Y (multiple visits are possible). If the set of the visited vertices is V C V , then we 
denote the family by Y pa? /,(V), and we say that V is the vertex set of the family. In a similar manner 
we can define the family of directed cycles, and denote it by T cyc i e (V), and say that V is the vertex set 
of the family. 

Definition 4.3. Let / : % C R" ->• W, Q) C %, and 5? be a cover of 0. We say that the directed 
graph <£(V,$) is a graph representation of f on 3> with respect to y, if there is a i : V — > 5? 
bijection such that the following implication is true for all H,v£f: 

/(i( M )n^)ni(v)n^0=> (u,v) e#, 

and we denote it by oc (/ ; ^ 5 

Having a graph representation ^ of / on @ with respect to , we take the liberty to handle 
the elements of the cover as vertices and vica versa, omitting the usage of i. It is important to 
emphasise that in general (u,v) G $ does not imply that f(uC\@) nvfl @ / 0. If we have (u,v) G 
# 44> f(u fl ^) fl v n *2l / 0, then we call an exacf graph representation. 

Enclosure of the non- wandering points 

Consider the continuous map 

/:%C1"^E" (4.1) 
Let f-\x) = {y£@ f : f(y) = x}, for x G W. 

Definition 4.4. The point q G ^/ is called & fixed point if /(#) = ^r. ^ G ^/ is a periodic point with 
minimal period m if f m (q) = q and for all < k < m : f k (q) / q; q G S>f is eventually periodic if it is 
not periodic, but there is a &o such that /** (q) is periodic. The point q G 3>f is a non-wandering point 
of / if for every neighbourhood U of q and for any M > 0, there exists an integer m> M such that 

/'«(£/ n%)n£/n%/0. 

Let A" C ^ be a compact set. We denote the set of periodic points of / in K by Per (/; K) and the 
set of non-wandering points of / in K by NonW (f;K). 



17 



Instead of directly studying the map ( |4.1[ ), we will analyse different graph representations of /. 
Let K C be a compact set satisfying /(Z) C ^, ^ a partition of K, and ^ oc (f,K, &>). We shall 
use the algorithm from [9] to enclose the non-wandering points in K. 

Algorithm 1 Enclosure of non-wandering points 
1: function ENCLOSURENONW(/,A', Sq) > /is the function, K is the starting region. 



2: fc<-0 

3: >6 «- Partition^, Sq) o ^) is a partition of K, diam(^) < So. 

4: loop 

5: +- Transitions(%,/) > The possible transitions (extra edges may occur). 

6: % <- GRAPH(^ ,4) >^«(/,d,^) 

7: for all v e ^ do 

8: if v is not in a directed cycle then 

9: remove v from §^ 

10: end if 

ll: end for 

12: if is empty then 

13: return > There is no non-wandering point in K NonW (/; A") = 0. 

14: else 

15: 4+1 «- 4/2 

16: %+i «- Partitioned, 5*+i) > is a partition of d|, diam(^ + i) < 4+1 

17: k^k + l 

18: end if 

19: end loop 



20: end function 



However, this algorithm for enclosing non-wandering points appeared without a full proof. We 
will give it here, not just for the sake of completeness, but because some steps are non-trivial. We 
need to take special care when a non-wandering point is on the boundary of a partition element. 

For any x E K, define 

# x :={ae^:j£«}. 
Since we are working with finite covers, both & and 3P X are finite. 

Lemma 4.5. For every q E K, there is an r\ q > such that for any u E if uC\B(q\r\ q ) ^ 0, then 
q Em. 

Proof. Since we are working with finite partitions, this is very easy to see. If & q = & then any 
positive number satisfies the condition. Otherwise define 

r/ := min d(q,u). 

This is positive, since & \ @P q is a finite set and for \ 3P q , u is compact and q £ u. Now we 

can take any number for r\ q from (0, tj ). □ 
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Lemma 4.6. For every q G NonW (f;K), there are u, v G ^ swc/z 

1. ^£«flv 

2. /or every e > 0, f/zere are x = x(e),N = N(x(e),e) G N toc/j f/iaf lim e ^oN(-*(£),£) = °°, * G 
B(q;e),f N (x) G B(q;s), x G w and f N {x) G v. 

Proo/ Consider a decreasing sequence of positive numbers {£/t}~ =0 with lim^oo £/t = and Eq < r\ q . 
Since o is non-wandering, 

3N k >k: f N *{B{q;e k ) n %)nB( ?;£t ) n % / 0. 

Therefore 

3x* G B(#fifc) : f k {x k ) G B (>;£*). 

Since ^ is finite, we can pick u G 2P q such that it contains infinitely many x^-s. We may assume, by 
switching to an appropriate subsequence and reindexing, that \fk : G u. It is now possible to choose 
- because of finiteness again - an element of the partition v G 3? q such that it contains infinitely many 
of f Nk (xk). Switching again to the subsequence, the required conditions are now satisfied. □ 

Remark 4.7. If there exists such that q G int(w) then it follows from the definition, that u and 

v = u is a good choice. By int(w) we mean the interior of the set u. 

Lemma 4.8. For every q G NonW (/;.£) \Per (f;K) there is an element u G & q and a family of 
directed cycles T cyc i e (V) in such that u G V, and the family encloses infinitely many trajectories in 
K of the form (x k ,f(x k ),f 2 (x k ),...,f Nk (x k )) with 

lim N k = °° and lim x k = lim / * (x k ) = q. 

k— >°° k— >°o k— >°° 

Proof. If there is a w G such that o G int(w), then since q G NonW (f',K), there are infinitely many 
such trajectories, and each one of them is enclosed by a directed cycle that passes through u. Since 
there is only finite number of families of directed cycles, therefore we can pick one family T cyc i e (V) 
that encloses infinitely many trajectories, and u G V. 

If such u cannot be found, then we will do the following: for a series of positive numbers 
liim^oo e k = 0, with the use of Lemma 4.6 we obtain the sets u,v G P q , the points x k G u,f Nk (x k ) G v 



with N k — > oo, and that lim^ooXjt = lim^oo f Nk (x k ) = q. 

Since {f Nk ~ Y (x k )} C K, and K is compact, we may assume that this sequence converges to a point 
q' G K, that is f Nk ~ l (x k ) — > q'. From the continuity of / it follows that 

f(q') = hmf Nk (x k )=q. 



Since & q i is finite, and because of Lemma 4.5 infinitely many of the points f Nk 1 (x k ) are inside one 



of its elements, without loss of generality we may assume, that 

VkeN:f Nk -\x k ) G v G & q <. 
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With similar argument we may assume that 

WeN:/(x t )e«'e% 

We are thus considering directed paths in the graph of the following form: 

u — y ii — \ . . . — > v' — y v. 

Since there is only finite number of families of directed paths in W and infinitely many trajectories 
of the desired property are enclosed by them, there must be at least one such family Tpat^iV) that 
encloses infinitely many trajectories itself. The directed paths V — > v — > u' and V — > u — > u' are 
present in , since f(q') = q. The situation is depicted in Figure|4] 



/ \ 




u v 



Figure 4: The directed paths close to q 

We can conclude our argument that the set V is a vertex set for a family of directed cycles as well, 
since 

U — Y U — y . . vertices fiomV . . — >• v' — >■ V — y U — y . . . vertices from 17 . , . v > M) 

is a directed cycle. Thus T cyc [ e (V) is well defined and encloses infinitely many trajectories of the 
desired type. □ 

Now we are ready to prove the correctness of Algorithm [T] 

Theorem 4.9. If NonW (/; K) ^ 0, then Algorithm^will never stop and NonW (f;K)Ci / k is satisfied 
for every k. 

Proof. Assume that q € NonW (f;K). If q is a periodic point with period N + 1 , then its orbit 
{<?, /(<?), • • • ,f N (q)} is enclosed in a directed cycle and none of these vertices is removed in the first 
step. Therefore q is a periodic point of / restricted to | f{ \ as well. Repeating the same argument gives 
that q is always enclosed. 

If q £ NonW(/;^), but it is not periodic, then we obtain a family of directed cycles T cyc i e (V) 
from Lemma pT8j that encloses infinitely many trajectories of the type mentioned before. The vertices 
in V are not removed, since they are in a directed cycle, thus the enclosures for all of these trajectories 
are preserved. Therefore q is a non- wandering point of / restricted to |^i| as well. Since |%| is 
compact, we can repeat the same argument and obtain families of directed cycles T cyc i e (Vk), that 
enclose infinitely many of these trajectories, that ensure that q stays a non-wandering point when we 
restrict / to |%| and q G |%| is satisfied. □ 
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Remark 4.10. The theorem does not imply that at one step each vertex containing q is kept. 

Remark 4.11. We decide whether a vertex is in a directed cycle by decomposing the graph into 
strongly connected components. The vertices that form a component by themselves and have no 
self edges are the ones that are not in directed cycles. To find this decomposition, we will use the 
algorithm of Tarjan l23l . that runs in linear time. 

Inner enclosure of the basin of attraction 

Consider the continuous map 

/:%C1"^1". (4.2) 

Definition 4.12. A set G C <&f is called invariant if f(G) = G. An invariant set G is called an 
attracting set if there exists an open neighbourhood % C <2)f of G such that 

(V open neighbourhood V D G) (3L(V) > 0) such that \fk > L(V) : f{%) C V (4.3) 

This neighbourhood *W is called a fundamental neighbourhood of G. The basin of attraction of G is 

UteN r*(^). 

Assume now that K C ^ is a compact set such that C X". When analysing the forward 
orbits starting from K we can work with f\g. Assume that there is an attracting set G C K for 
and a neighbourhood £/ such that G C U C K and C7 is inside the basin of attraction of G. We want 
to find a set B in the basin of attraction of G so that U C B. We will use the algorithm from : 

Algorithm 2 Inner enclosure of the basin of attraction 
l: function Basin_of_Attraction(/,^, So; U) > C/ is the initial attracting neighbourhood. 



2: fc<-0 

3: W > We collect the vertices in the basin of attraction into W. 

4: % <- Partition^, do) t> ^) is a partition of K, diam(^o) < 8q. 

5: loop 

6: Sk ^— Transitions(^ UW,f) > The possible transitions (extra edges may occur). 

7: <S k <r- GRAPHS UW, 4) >9k°<(f,\%UW\,Y k \JW) 

8: repeat 
9: ready i- TRUE 

10: for all v G ^ do 

ll: ifvC£/u|W|or/(v)C£/u|W|then 

12: move v from % to W o v is contained in the basin of attraction of G 

13: ready <- FALSE 

14: end if 

15: end for 

16: until ready > The remaining vertices are not attracted at this resolution 
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Algorithm 2 Inner enclosure of the basin of attraction (continued) 



17: if STOP(fc, %,W, 4) then > Some stopping condition 

18: return W 

19: end if 

20: 4+1 <- 4/2 

21: %+\ <- Partition(|^|,4+i) > %+i is apartition of \%\, diam(%+i) < 4+1 

22: it f- it + 1 
23: end loop 



24: end function 



In Algorithm|2| we move those partition elements into W, that lie inside or are mapped into the 
initial attracting neighbourhood, or the other elements in W. Then we refine our remaining partition 
and continue our procedure with a new one, that has a diameter half as big as before. Since at the 
beginning W was empty, it will only contain sets that are contained in the basin of attraction of G. 
Thus after each cycle, |W| is an inner enclosure of the basin of attraction of G. We stop our iteration 
when a stopping condition is satisfied, for example 4 < A, where A is a small positive number given 
in advance. 



5 The completion of the proof of Theorem 1.1 



Consider now a parameter value a for F a such that 0.5 < a < 1. In order to prove that the fixed 
point a is globally attracting, we need the following observation: given any starting point (xo,yo), 
the accumulation points of the sequence {F^(xo,yo))T_ are non- wandering points of F a . We want to 

show that the only non-wandering point of F a in is the fixed point. We know from Lemma 2. 1 

(a) 

that it is enough to look for points in S] ,i & No. Thus our goal is to prove that 

la) 

1. Sj is entirely in the basin of attraction of a, 

(a) 

2. or equivalently, S] contains exactly one non-wandering point, and that is a. 

Our strategy is to divide the parameter range [0.5, 1] = [0.5,0.875] U [0.875,0.999] U [0.999, 1] into 
small subintervals [a]. The diameter of these subintervals will vary between 10~ 3 , 10~ 4 and 10~ 5 in 
practice. For one small parameter interval [a] we shall follow these steps: 

1. Let i > 1 be the smallest integer such that \hf® ]) -^■ ( J + 1 1 ) |+|gl ( | C<1) l< 10_9 - 

2. The function ConstructRegionQa]) returns a rigorous enclosure [S] such that 

U ^ C [5]. 



3. Using Propositions 3.1| 3.2 and 3.3 the function FindAttractionDomain([a]) returns an 



£o > such that K(a; £q) is contained in the basin of attraction of a for every a £ [a]. 
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4. Enclose rigorously U ae [ a ]NonW (F K ; [S]) \ {«}) by removing parts of [S] that do not contain 
non-wandering points of F a or are in the basin of attraction of the fixed point a for every 
a £ [a]. We do this by simultaneously checking the criteria from line[8]of Algorithm [I] and line 
1 1 of Algorithm[2| If we obtain an empty enclosure at some step, then we have proved that the 
fixed point in the given parameter region is globally attracting. 

We sum it in the following algorithm: 



Algorithm 3 Proving the global stability of a for the Ricker-map 



l: procedure RiCKER([a], 5) 

2: [S] <— ConstructRegion([a]) 

3: £o FindAttractionDomain([a]) 

4: [U] <-K([a];eb-(a + - or)) 

5: r<- Partition([5], 8) 

6: repeat 

7: g<- Transitions^, F [a] ) 

8: £f ^GRAPH(r,<r> 

9: T <— {v : v is in a directed cycle } 

10: for all v e f do 

11: if v ^ T or v C [f/] or F [o] (v) C [17] then 

12: remove v from ^ 

13: end if 

14: end for 

15: 8^8/2 

16: Partitioned, 5) 

17: until \y\=% 
18: end procedure 



> from (2.1 1 



> from Propositions 3.1 3.2 and 3.3 



> V is a partition of [5], diam^) < 8 
> The possible transitions (extra edges may occur). 

>sfo e (F [0] ,m^) 

> with the use of Tarjan's algorithm 



We know that FindAttractionDomain[a] returns an £o > suc h that for every a G [a] 
[a _ ,a + ], the set K(a;£o) is in the basin of attraction of a. Assume that a + — a - < £q and let 



; £o - (a + - a 



Now £ > and the set K([a];e) is in the basin of attraction of a for every a G [a]. Observe that, 



using subintervals with or — a < 10 -3 and the £o obtained from Propositions 3.1 3.2 and 3.3 
or — or < £o is satisfied. 

After each step in the main cycle in Algorithm [3j \Y\ is a rigorous enclosure of all the non- 
wandering points of F a in [5] \ {a} for every a E [a]. This is easy to see since vertices are removed 
for two possible reasons which are both checked in line 1 1 of Algorithm [5] First if v ^ T, then v does 
not contain non- wandering points for any F a , a € [a] as we have seen in the proof of the correctness 
of Algorithm[T] Second if v C [U] or F [a] (v) C [U] that is v is inside or mapped into the small attracting 
neighbourhood of every fixed point. Note that if a vertex is inside the basin of attraction of a fixed 
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point a, then it cannot contain any other non- wandering point of F a , not even on the boundary. This 
is a similar criterion to what we have used in line 1 1 of Algorithm [2] The difference is that now we 
remove these vertices, consequently we do not have to collect them into a list. If the procedure ends 
in finite time, then we have established, that there are no other non-wandering points in [5] , thus the 
fixed point is globally attracting for all a G [a] . We state this as 

Proposition 5.1. If Algorithm^ends infinite time with input parameters ([<x],10 _1 ), that is, after 
finite number of steps, \¥\= is satisfied, then a is a globally attracting fixed point of the two 
dimensional Ricker-map F a for every a £ [a]. 

We implemented our program in C++, using the CAPD Library [51 for rigorous computations, 
and the Boost Graph Library [22] for handling the directed graphs. The recursion number in Tarjan's 
algorithm was very high, therefore we converted it into a sequential program, using virtual stack 
structures from the Standard Library in order to avoid overflows. Instead of simulating the Ricker- 
map itself, we used its third iterate, the formula is still compact enough not to cause big overestimation 
in interval arithmetics and it considerably speeds up the calculations. 

As an example, we ran our program for the parameter slice [0.9,0.90001], with 8 = 10 > 1 as the 
initial diameter for the partition. We show the evolution of the enclosure during the first 8 iterations 
on Figure [5] 








Figure 5: Enclosure after 1, 3, 4, 6, 8 steps 



The small rectangle is the attracting neighbourhood [U ] . After 6 iterations the diameter of the partition 
is sufficiently small in order to have some boxes removed from the inside even though they are in 
directed cycles. This happens because they are contained in, or get mapped into the basin of attraction 
of the fixed point. 

We used different sizes for the parameter intervals and ran the computations on a cluster of the 
NIIF HPC centre at the University of Szeged (48 cores, 128 GB memory / cluster) parallelising it 
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with OpenMP. We summarise some technical details in Table [T] 



parameter 


size of slices 


# of CPU 


max. memory 


wall clock time 


total time 


[0.5,0.875] 


10" 3 


48 


2.30 GB 


1.2s 


38.5s 


[0.875,0.95] 


10~ 3 


48 


3.05 GB 


2.4s 


52.3s 


[0.95,0.99] 


10" 4 


48 


3.07 GB 


33.7s 


22m 5.9s 


[0.99,1] 


10" 5 


20 


65.30 GB 


204m 42.6s 


1800m 37.1s 



Table 1 : Resources used by the program 



Remark 5.2. Here wall clock time (real) refers to the actual running time of the process, whilst the 
total time (user + sys) is the sum of the time spent on individual CPUs. 

The complexity of the computations for some parameter slices is shown in Table [2] 



parameter slice 


[S] 


# iteration 


max # of vertices 


max # of edges 


[0.875,0.876] 


[2.072e-04,5.049031] 2 


13 


242 


1,676 


[0.999,0.99901] 


[2.928e-06,7.369087] 2 


27 


729,528 


4,193,329 


[0.99999, 1] 


[2.822e-06,7.389015] 2 


33 


3,105,304 


118,751,916 



Table 2: Complexity of the computations 



The program ran successfully, thus we established that the fixed point is globally attracting for 
(X £ [0.5, 1]. Combining this with Proposition 2.2 and Proposition 3.4 the proof of Theorem |1.1| 
is completed. 
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A Appendix 



, , 1 / . , — — \ / / , z-izs/-l+4a+z+iy/-l+4az \ , , 

s(^a) = —==(- l + ^l+4a)((-l + e * Uz + z) 



(A.l) 

- izV-l+4a+i+iV-l+4a; 2z 2lZ \ \ 

+ a[-l + e 2« -| -| 

1 + «V-1 +4a i + V-l +4a / / 



^ o(a) ^4 + ^m« (A - 2) 

v ' -i+4o + iV-l+4o 

£02 a = V . . ' (A.4) 

v 1 V- 1 + 4a (i - 2ia + V- 1 + 4a) 

, , i(l + a-iV-l+4a) 
av-l +4a 

, . -3/ + 2/a + V-l+4a 

£21 (a) = — . . : . (A.6) 

2aV-l +4a 

, , 2(3;'-2;a + V-l+4a) 
8n(a) = / - V 7 _ ; 2 (A.7) 
V-l+4a (/ + v/-l+4a) 

, . 4;'a (5/ + V-1 +4a) 
£03 (a) - - 1 / - \ 3 ( A - 8 ) 
V-l+4a (i + V-l+4a) 

, N 8(-2 + a-2?V-l+4a) 
£4o(a) = - 1 \ 3 (A.9) 

, . 2(-2/ + /a + V-l+4a) 
g3l(a) = -a H + 4^aW-l + 4a) (A ' 10) 

g22 (a) = ? 2{ - 2 + a l , (All) 

' o(-l+4o + iV-l+4a) 

, . 8 (2- a- W-l+4a) 
813(a) = — ; ^ (A. 12) 

V-l+4a (/ + V-l+4a) 

£ 04 (a) = / 7 - ,4 (A.13) 

V-l+4a (? + V-l+4a) 
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hnv.a 



(z) = z- /i20 ^ + I z 3 (3/ !20 (a) 2 -/130(a) + 3/ZH (a)/z 02 (a)) + ^ 4 (-15%(a) 3 + 10%>(a)%)(a) 
- 30/jh (a)feo(a)^02(a) - 3/z 2(a)Ao2(a) 2 + 4/jn(a)A 3(a) - 12An(a)ft 2(a)Aii(a)) 

- An (a)zz+ ^ 2 (3An(a)A 2 o(a) + A 02 (a)Ao 2 (a) + 2Aii(a)An(a))z 

+ ^z 3 (-15/!ii(a)/!2o(«) 2 + 4/z 11 (a)/z 3 o(a) - 12h n {a) 2 h 02 (a) + 3h l2 (a)h 02 {a) 

- 6h 02 (a)h 2Q (a)h 02 {a) + h Q2 {a)h 03 (a) - 12h n (a)h 20 (a)h n (a) - 6h 02 (a)h 02 (a)h n (a) 

- 6h n (a)h n (a) 2 + 3h n {a)hi 2 (a) - 3/jn(a)/j 2(a)/j 20 (a)) z 

- ho2 ^ z + ^z(2h n (a) 2 - h n {a) + h Q2 {a)h 2 Q(a) + 2h 02 (a)h u (a) + h u (a)h 20 (a)^z 2 

+ ^z 2 ^-12/!ii(a) 2 /z 20 (a) + 3/zi 2 (a)/! 20 (a) - 3/z 02 (a)/i 20 (a) 2 + A 02 (a)A 30 (a) + /! 03 (a)/zo2(a) 
- 9ho 2 {a)hn(a)ho 2 (a) - 12/zn(a) 2 /zn(a) +4/n 2 (a)/m(a) - 6h 2(a)h 2 o(a)h n {a) 



- 6h 02 (a)h n (a) + 2A 02 (a)/i 12 (a) - 3Aii(a)A 20 (a)A 20 (a) - 3/z 02 (a)/z 02 (a)/i 20 (a) 
- 6/!ii(a)/!n(a)/i 20 (a))z 2 + l - (-h (B (a) + 3h Q2 (a)h n (a) + 3h 02 (a)h 2Q (a)^ z 3 

+ ^z(-6/zn(a) 3 + 6h n (a)h l2 (a) + h Q3 (a)h 20 (a) - 9A 02 (a)An(a)A 20 (a) - 3/z 02 (a) 2 /z 02 (a) 

+ 3/z 3(a)/in(a) - lSho 2 (a)hn{cc)h n {cc) - 6h n (a) 2 h 2 o(a) + 3Ai 2 (a)A 20 (a) 
- 3/!o2(a)/! 2 o(a)/z 20 (a) - 12Ao 2 (a)Aii(a)A 20 (a) - 3/!n(a)/z 20 (a) 2 + /zn(a)/i 30 (a))z 3 

+ (4/! 03 (a)/zii(a) - 12/! 02 (a)/zn(a) 2 + 6/! 2(a)Ai2(«) - 3/i 02 (a) 2 /z 20 (a) - 12/z 02 (a) 2 /zn(a) 

+ 6/zo3(a)/i 20 (a) - lSho 2 {a)hn{cc)h 20 (a) - 15h 02 (a)h 2Q (a) 2 + 4/zo 2 (a)/z 3 o(a))z 4 

(A. 14) 

, 2(l + x/l -4a-a) 

A20 a = — - ; , (A.15) 

y ' o(l-4o + iV-l+4a) 

*n(a) = - 7 2( ~ 1 + a) T (A.16) 

v ; o(-l+4o-iV-l+4a) 

, N 2a (3/ + V-l+4a) 

hoi(a) = 2 ^ ' ; (A.17) 

(i + V-I +4a) (-i + 4/a + aV-l +4a) 

/130(a) = 

(l6 (-12 - 12/V-l +4a + a (45 + 2W-1 +4a + a f-27 + 2a - 5iV-l +4a)))) 
/ 4 \-i (A-18) 

( (l + Vl - 4a) V-l +4a (i + V-l +4a + a (-2i + a (-5i + V-l + 4a ))) ) 

Ai 2 (a) = 

((4/ + a (-4 (5i + 3 V-l + 4a) + a (281 + 18V-1 +4a + a (-31* - 9V-1 +4a + 
a (9/ + V-l +4a)))))) ((a 3 (-i + V-l + 4a + a (6/ + a (-13/ - 15V-1 + 4a + (A. 19) 
2a(lO/ + V-l+4a))))))~ 1 

h m {a) = 
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( (2 (5i + 7V-1 + 4a + a (-29/ - 14V-1 + 4a + a (l7i - ia + 3V-1 + 4a) ) ) ) ) 
( (/ + V-l + 4a + a (-2 (7i + 6\/-l + 4a) + a (70/ + 46V- 1 + 4a + 
a (-6 (23/ + 10V-1 + 4a) + a (73/ - 4/a + 13V- 1 + 4a)))))) ' 



|^(a)| + |^ v (a)| + |C(«)| 



-2(-l + a)a + ^(2 + a) + ^ ,^(4 + a) J (V ^" 1 +4o >) 
|C(«)| + |C(«)| + |^(a)| + |C(«)| = 



1 36+ 15a- 12a 2 + 4a 3 /6 - 40a + 75a 2 - 21a 4 + 4a 5 + 4a 6 
6~y a 5 (-2 + 7a + 4a 2 ) + y (2 - 15a + 22a 2 + 23a 3 + 4a 4 )4a 6 

1 / 12 - 54a + 66a 2 - 13a 3 + 4a 4 + 4a 5 
+ 6 y a 3 (-1 + 13a - 57a 2 + 88a 3 - 18a 4 + 7a 5 + 4a 6 ) 

1 h - 16a + 13a 2 + 124a 3 - 198a 4 + 15a 5 + 54a 6 + 9a 7 
2y a 6 (-l +4a)(-l +4a + a 2 ) 2 

lC(«)l= 



^- f — 1140 + 10874a - 12660a 2 - 143073a 3 + 423177a 4 
24 V 

- 211261a 5 + 1356a 6 + 30869a 7 - 13766a 8 - 2238a 9 + 1489* 
+ 256a 11 ) 2 ((1 - 4a) 2 a 6 (-1 + 4a + a 2 ) 3 (-1 + 5a - 2a 2 + 



a 



i 

3\\ Z 



lC(«)l= 

(-36 + 912a - 9354a 2 + 49026a 3 - 133548a 4 + 155248a 5 + 24851a 6 - 182342a 7 
+ 127014a 8 - 47122a 9 - 12543a 10 + 34528a 11 - 329a 12 + 1925a 13 + 2452a 14 
+ 361a 15 ) 2 (-l+5a-2a 2 + a 3 )~^ (6a 5 (2 - 15a + 22a 2 + 23a 3 + 4a 4 )) ' 

lC(«)l= 

(- 16 + 4a + 2240a 2 - 14868a 3 + 19782a 4 + 91297a 5 - 309731a 6 + 259610a 7 
+ 81080a 8 - 147591a 9 + 12815a 10 + 19871a 11 - 8236a 12 + 241a 13 + 2132a 14 
+ 361a 15 ) (-1 + 5a - 2a 2 + a 3 )" 1 (4a 5 (2 - 15a + 22a 2 + 23a 3 + 4a 



\hZ(a)\= 

^ ((36 - 840a + 7728a 2 - 35454a 3 + 84157a 4 - 96139a 5 + 39017a 6 
+ 15361a 7 - 22836a 8 + 10489a 9 + 5142a 10 - 397a 11 + 922a 12 + 529a 13 
+ 64a 14 )) 2 (a 9 (-l +4a + a 2 ) 3 (2 - 17a + 35a 2 + 4a 3 - a 4 + 4a 5 )) 
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IC(«)I= 

^- (36 - 396a + 1350a 2 - 1422a 3 + 318a 4 + 441a 5 - 145a 6 + 100a 7 

24 V , (A.27) 

+ 25a 8 ) 2 (a 6 (-l+4a + a 2 ) 2 (l-9a + 22a 2 -9a 3 +4a 4 )) 2 

k 2 4 V)l= 

^ (-432 + 2808a - 5016a 2 + 4692a 3 - 2584a 4 
24 V 

+ 930a 5 - 210a 6 + 7a 7 + 27a 8 - 6a 9 + a 10 ) 2 ( A28 ) 
(a 4 (-l +4a + a 2 ) 2 (l -9a + 22a 2 -9a 3 +4a 4 )) 2 

Vl\a)\= 

I (90 - 810a + 1635a 2 + 1485a 3 - 3462a 4 + 4056a 5 - 659a 6 
6 V 

- 108a 7 - 320a 8 + 24a 9 + 87a 10 - 102a 11 + 36a 12 + a 13 ) 2 ( A29 ) 
(a 7 (2- 17a + 33a 2 + 13a 3 + 4a 5 + a 6 ) (-1 +4a) 2 ) 2 

Vf{a)\= 

1 f 64 - 608a + 1348a 2 + 728a 3 - 1692a 4 + 1948a 5 - 688a 6 - 2008a 7 + 1184a 8 
4 V 

+ 742a 9 - 355a 10 - 136a 11 + 34a 12 + 14a 13 + a 14 ) 1 ( A - 30 ) 
(a 9 (-l +4a)(-2 + 7a + 6a 2 + a 3 ) 2 ) 2 

l - (-54 + 918a - 5973a 2 + 18921a 3 - 32250a 4 + 32742a 5 - 15643a 6 + 2506a 7 

+ 3246a 8 - 3551a 9 + 1327a 10 - 156a 11 - 48a 12 + 59a 13 + 16a 14 + a 15 ) 2 (A31) 
(a 7 (-2 + 9a + a 2 + a 4 ) (6 - 48a + 90a 2 + 4a 3 ) 2 ) 

|r 2 4 V)l= 



^ (-432 + 2808a - 5016a 2 + 4692a 3 - 2584a 4 
24 V 



+ 930a 5 - 210a 6 + 7a 7 + 27a 8 - 6a 9 + a 10 j 2 ( A32 > 

(a 4 (-l +4a + a 2 ) 2 (l -9a + 22a 2 -9a 3 +4a 4 )) 2 

ci(a)= (A.33) 
(li - 2V-1 +4a + a (2 (-7i + 5V-1 +4a) 



+ a (25i - 13V-1 +4a + a (-25i + 7V-1 +4a - a (-7? + V-l + 4a ))))) 
(2aV-l + 4a(-/ + \/-l + 4a) (/a + V-l + 4a)) 
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